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Представлена вычислительная технология, обеспечивающая доступ к информационным материалам в области дендроэкологии 
и программным средствам для анализа изображений годичных колец деревьев с целью выявления изменений параметров окру- 
жающей среды, отраженных в приросте дерева. В основу моделирования годичных колец деревьев положены фундаменталь- 
ные свойства радиального роста дерева - монотонность роста во времени и пространстве и ограниченность его скорости, а из- 
менения плотности древесины рассматриваются как пространственно-временной колебательный процесс. Технология базиру- 
ется на динамическом управлении распределёнными вычислительными ресурсами. 


1. Введение 

При реализации достаточно сложных задач 
очень быстро проявляется недостаток имеющихся 
вычислительных мощностей, заключающийся, как 
правило, в увеличении времени отклика системы. В 
частности, особенность задач дендроэкологическо- 
го анализа заключается в необходимости математи- 
ческой обработки очень большого объема данных 
(иногда требуется обработка временных рядов на- 
блюдений за 1000 лет и более), широким спектром 
решаемых задач, коллективной (одновременной) 
работой многих сотрудников, практически, на всех 
этапах дендроэкологических исследований, а также 
в необходимости хранения и систематизации боль- 
ших объемов неоднородной структурированной 
информации (собственно хронологические ряды 
наблюдений, результаты обработки, сопутствую- 
щие метеорологические, геологические, геофизи- 
ческие, аэрокосмические и т.п. ряды наблюдений). 

Отмеченные выше особенности задач дендроэ- 
кологии не позволяют создавать профессиональ- 
ные центры обработки и хранения дендроэкологи- 
ческих данных в отдельно взятом научном учрежде- 
нии. Слишком большие вычислительные возмож- 
ности и мощные базы данных требуются для этого и 
слишком велики материальные затраты. Естествен- 
ным решением в сложившейся ситуации является 
распределение информационного, программного и 
технического обеспечения (ресурсов) решения за- 
дач дендроэкологии в нескольких специализиро- 
ванных центрах обработки и виртуализация доступа 
к этим информационно-обрабатывающим ресур- 
сам исследователей. Такая концепция однозначно 
подразумевает в качестве одного из системных ре- 
шений использование в таких профессиональных 
центрах технологии распределенной обработки 
(^«/-технологии), обеспечивающей динамическое 
изменение основных компонентов инфраструкту- 
ры системы обработки (от структур хранимых дан- 
ных - до схем и алгоритмов решаемых задач) и по- 
вышение вычислительной мощности. 

Наиболее значимым при решении задач ден- 
дроэкологии становится восстановление ретро- 
спективных сведений об изменчивости условий 


окружающей среды, при изучении которых ис- 
пользуются результаты, полученные путём инстру- 
ментальных измерений. Такая информация являет- 
ся наиболее ценной и широко используемой, но 
часто обладает существенными недостатками, ко- 
торые ограничивают её использование, в частно- 
сти, кратковременность и нерегулярность прово- 
димых наблюдений, их территориальная ограни- 
ченность и т.д. В связи с этим перспективным под- 
ходом к реконструкции хода природных процессов 
являются индикационные методы, основанные на 
наличии связей между составом и структурой при- 
родных тел и параметрами внешней среды. Инди- 
каторами в таких методах обычно являются объек- 
ты живой природы, характеризующиеся неодно- 
родными анатомо-морфологическими характери- 
стиками, возникшими в разные годы, сезоны года 
и даже в течение суток, и сохраняющими такое 
строение в течение длительного времени. Наибо- 
лее известным и широко распространенным при- 
мером регистрирующих структур служат годичные 
слои (годичные кольца) древесины в стволах, ве- 
твях и корнях деревьев и кустарников [1-3]. Отме- 
тим, что годичные кольца дают интегральную ин- 
формацию о качестве жизни и в этом их преимуще- 
ство перед физическими методами измерения па- 
раметров окружающей среды, которые могут с вы- 
сокой точностью определить значения конкретных 
параметров (концентрацию веществ, интенсив- 
ность излучения, скорость движения, температуру 
и прочее). Однако живому организму важно их оп- 
тимальное сочетание, что и фиксируют годичные 
слои древесины [3,4]. 

Таким образом, исследование особенностей 
роста дерева (от начала до завершения) с деталь- 
ным изучением ростовой активности внутри сезо- 
на с использованием развитого математического 
описания изменения профиля плотности древеси- 
ны, поддержанное численными алгоритмами, реа- 
лизованными в высокопроизводительных вычи- 
слительных системах с новыми эффективными ра- 
спределёнными методами обработки дендроэколо- 
гических данных, даст возможность получить но- 
вые результаты и качественные оценки экологиче- 
ских параметров окружающей среды. 
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2. Сравнительный обзор и анализ 
аппаратно-программных средств обработки 
и анализа годичных колец деревьев 

Анализ научных исследований в области ден- 
дроэкологии, дендрохронологии и дендроклимато- 
логии показывает, что в настоящее время для моде- 
лирования и анализа дендроэкологических данных 
используется достаточно широкий спектр аппарат- 
но-программных средств [5, 6\. Использование 
программных средств зависит от задач исследова- 
ния. Это может быть статистическая обработка 
первичных данных замеров характеристик прирос- 
та и получение надёжных обобщённых хроноло- 
гий, сопоставление характеристик прироста с фак- 
торами внешней среды, моделирование процессов 
роста, изучение пространственного распределения 
характеристик прироста и визуализация результа- 
тов анализа. Для решения этих задач одни исследо- 
ватели предпочитают использовать широко рас- 
пространённые статистические пакеты (8А8, Маі- 
ІшЬ, 8у$1а1, 8ТАТІ8ТІСА и др.), электронные та- 
блицы (СШАТТКО, ШТШ, Ехсеі и др.), универ- 
сальные системы математической обработки ре- 
зультатов, предназначенные для численного и сим- 
вольного решения математических задач различ- 
ной сложности (МагЬСАЭ, МаШетаІіса и др.) и 
интегрированные программные решения такие как 
РКЕССЖ 5.17С, ТЯЕЕКШО 3.0, Еі§поѴшоп 
1.32, Т$АР-\Ѵіп Ргоіекыопаі 0.30, ОепсІгоСІіт 2002, 
\УіпОЕШЯО, ОепбгоЕаЬ 470, РА8Т 32, 08М 3.10, 
Меазиге12Х, или решения на базе географических 
информационных систем (Е8Я1 АЯС/ШЕО, Аг- 
сѴіе\ѵ, Маріи Го, АиІоСАО Мар и др.) с форматами 
различных баз данных (Місго$ой Ассеж. Огасіе, 
СІВА8Е, ЕохРго и др.) для связи обрабатываемых 
данных с географическим местоположением их 
сбора и отображением на электронных картах [3]. 

Рассмотренные программные комплексы и си- 
стемы достаточно универсальны, ориентированы 
на широкий класс исследований и выполнение ча- 
сто используемых операций обработки и анализа 
дендроэкологических данных (например, измере- 
ние ширины годичных колец и др.). Однако, они 
плохо учитывают специфику требований конкрет- 
ных исследователей (например, определение ази- 
мута и среднеквадратичной ширины области мак- 
симального прироста и др.). 

3. Монотонность и непрерывность роста как основа 
построения математической модели годичных колец 

Математическая модель радиального сечения 
дерева строится, исходя из формы и структуры 
концентрических годичных слоёв, которые форми- 
руются на поперечном спиле дерева в процессе его 
роста. Для построения модели вводится полярная 
система координат, начало которой совпадает с 
биологическим центром диска дерева. Годичные 
кольца рассматриваются как функция двух пере- 
менных: р (радиуса) и Ѳ (азимута). На рис. 1 пока- 
зано поперечное сечение ствола дерева (пихты си- 


бирской, АЫез хіЫгіса ЬеАеЬ), полученное в резуль- 
тате сканирования с необходимой дискретностью, 
и профиль изменения интенсивности окраски 
вдоль радиуса сканирования (функция Я{р), отме- 
ченной пунктиром). На графике значения функ- 
ции Я(р) представлены в относительных единицах, 
а аргументом является расстояние вдоль радиуса 
сканирования. Заметим, что фрагмент значений 
функции Я(р) представлен по наиболее наглядно- 
му участку радиуса сканирования. 




Принимаем функцию Я(р) близкой функции 
плотности древесины и будем называть её функци- 
ей плотности. Функция Я(р) представляет собой 
колебательный процесс, период которого (период 
радиального роста дерева) непостоянен. В конце 
каждого периода, на границе между поздней и ран- 
ней древесиной, происходят резкие изменения 
функции плотности, которыми и фиксируется вре- 
мя (г) при соответствующих значениях р, что фак- 
тически определяет на радиусе временные реперы. 
Вследствие того, что радиальный рост не может 
убывать и меняться скачком во времени в пределах 
вегетационного периода, считаем, что функция 
роста р(т) является монотонной и имеет ограни- 
ченную производную. По этой причине обратная 
ей функция ср(р) также монотонна и, более того, 
она является строго возрастающей функцией. 

Из экспериментальных данных следует, что се- 
чение Я(р) содержит шумы, наложенные на основ- 
ное колебание. Локальный импульсный шум И(р) 
определён в области ранней древесины. Этот шум 
вызван главным образом недостаточной чистотой 
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обработки поверхности поперечного спила и есте- 
ственными флуктуациями размера клеток древеси- 
ны. В виду того, что формирование клеток проис- 
ходит под действием большого числа различных 
факторов, данный компонент имеет случайный ха- 
рактер. Локальный импульсный шум занимает 
внутреннюю область годичного кольца, имеет ма- 
лую амплитуду, по отношению к амплитуде полез- 
ного сигнала, и расположен в области высоких ча- 
стот. Наряду с высокочастотным шумом, имеется 
глобальный низкочастотный шум І(р), связанный с 
наличием возрастных зон на поперечном спиле 
ствола дерева: сердцевины и заболони. Принима- 
ем, что этот компонент шума накладывается на по- 
лезный сигнал в виде колебания, период которого 
значительно больше периода, соответствующего 
наибольшему годичному кольцу. 

Исходя из вышесказанного математическую 
модель радиального сечения годичных колец мож- 
но представить в виде суммы шумов и действитель- 
ного колебания Щр), модулированного по ампли- 
туде и фазе: 

К(р)=С+Щр)+Н(р)+Кр)> 0 , (*) 

1Лр)=а(р)ссыр{р), 

ср(р)=у~т{р). 


где о(р) - амплитуда колебания, ср(р) - фаза коле- 
бания, т(р) - время роста, Т - период вегетации, 
С - постоянная. 


Функция Я(р) формируется биологическим 
объектом под влиянием внешней среды и имеет вид 
интерференционной картины. Поэтому естественно 
предположить, что функция плотности древесины 
является результатом суммарного взаимодействия 
“биологической” и “экологической” составляющих. 
Отметим, что подобная формализация радиального 
сечения диска дерева позволяет эффективно приме- 
нить её для анализа как внешних, так и внутренних 
особенностей роста дерева в вегетационный период. 
Основой для проведения прикладного анализа роста 
являются свойства компонентов, входя щ их в по- 
строенную математическую модель. Главной и наи- 
более информативной составляющей модели явля- 
ется фаза, так как она связана с функцией роста об- 
ратно пропорциональной зависимостью. 

Для восстановления роста из выражения (*) на- 
ходится фаза, а затем определяется обратная ей 
функция. После последовательных преобразова- 
ний функция роста будет иметь вид: 


р(т)=г-‘ 



Функция роста р(т) даёт время появления каж- 
дой клетки дерева на исследуемом радиусе в про- 
цессе вегетации, т.е. для всех значений р в интерва- 
ле между реперами. Полученное аналитическое 
выражение уже пригодно для создания алгоритмов 
автоматического вычисления годового прироста 
деревьев за вегетационный период. 


Дальнейшие уточнения модели радиального се- 
чения диска дерева связаны с различными метода- 
ми определения амплитуды и фазы. В частности, 
широкое распространение для вычисления ампли- 
туды и фазы получил аналитический сигнал, широ- 
ко используемый в теории колебаний и волн. Ана- 
литический сигнал является комплексным, его дей- 
ствительная часть есть исходная функция, а мнимая 
часть образуется применением преобразования 
Гильберта к действительной части. Аналитический 
сигнал определяет способ измерения амплитуды, 
фазы и угловой частоты сигнала. Если спектр сиг- 
нала сосредоточен в относительно узком интервале 
частот (квазимонохроматический сигнал), то изме- 
нения амплитуды и фазы за время, соответствую- 
щее периоду основной частоты, при определенных 
условиях оказываются минимальными [7]. 

4. Структура и программная реализация 

многомашинного комплекса обработки 

дендроэкологических данных 

Авторами были проведены экспериментальные 
исследования требуемой вычислительной мощно- 
сти на такой задаче в области дендроэкологии как 
обработка изображений годичных колец деревьев 
(сканирование, обработка фильтрами). В зависи- 
мости от размера обрабатываемого изображения, 
требовалось до нескольких дес. мин процессорного 
времени (Сеіегоп 1,7 МГц, БОК ВАМ 256 Мб). На- 
грузка на процессор возрастала, практически, про- 
порционально числу подключившихся пользовате- 
лей. Одним из решений проблемы недостаточно- 
сти вычислительных ресурсов является кластериза- 
ция нескольких компьютеров в многомашинный 
комплекс обработки (рис. 2), так как реально уве- 
личивать вычислительную мощность одного ком- 
пьютера до бесконечности проблематично. 
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Рис. 2. Архитектура комплекса обработки дендроэкологиче- 
ских данных 
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Планирование и диспетчирование процессами 
обработки дендроэкологических данных возложе- 
но на специальный коммуникационный сервер. 
Взаимодействие всех пользователей с комплексом 
обработки осуществляется только через коммуни- 
кационный сервер. Непосредственно на этом сер- 
вере математические вычисления не выполняются. 
Основной задачей коммуникационного сервера 
является обеспечение оптимальной загрузки 
имеющихся в его распоряжении вычислительных 
серверов и обеспечение режима работы в реальном 
времени ( опііпе ) пользователей. Решение задач ма- 
тематической обработки данных (например, таких 
как вычисление радиального роста дерева и плот- 
ности древесины в период вегетации, проведение 
спорово-пыльцевого анализа последовательных 
слоёв торфа или построение длительных древесно- 
кольцевых хронологий и т.д.) возложено на вычи- 
слительные серверы кластера. 

Для информационной поддержки функциониро- 
вания многомашинного комплекса обработки ден- 
дроэкологических данных использована СУБД Огас- 
1е. С сервером базы данных Огасіе взаимодействуют 
как вычислительные серверы, так и коммуника- 
ционный сервер. Вычислительные серверы класте- 
ра, реализующие основную математическую обра- 
ботку, и коммуникационный сервер выполнены на 
базе АрасНе НТТР Бегѵег, а клиентская часть, реали- 
зующая интерфейс удалённого взаимодействия с 
пользователем, использует М§ Іпіегпеі Ехріогег. 

Рассмотренная выше схема кластеризации ре- 
сурсов была реализована для решения задачи ана- 
лиза изображений годичных колец деревьев. Вычи- 
слительный кластер обеспечивал решение таких 
задач, как сбор, систематизацию и хранение науч- 
ной информации в области дендроэкологических 
исследований, формирование и ведение базы экс- 
периментальных дендроэкологических данных, 
вычисление азимута и среднеквадратичной шири- 
ны области максимального прироста годичных ко- 
лец деревьев, вычисление значений индексов при- 
роста ширины годичных колец деревьев, оценку 
экологического состояния окружающей среды. 
Для указанных задач были разработаны и реализо- 
ваны авторские вычислительные алгоритмы. В ка- 
честве критерия, определяющего качество алгорит- 
мов и их программной реализации, было выбрано 
время выполнения. Оптимизация осуществлялась 
по минимуму времени выполнения [8]. 

Ниже приведена иллюстрация функционирова- 
ния комплекса при решении некоторых из перечи- 
сленных выше задач. В задаче вычисления азимута и 
среднеквадратичной ширины области максималь- 
ного прироста годичных колец деревьев исследова- 
ние включало: определение центра колец, сканиро- 
вание изображения по всем азимутам, преобразова- 
ние изображения из полярной системы координат в 
декартовую, обработку изображения медианным 
(рис. 3) и полиномиальным (рис. 4) фильтрами, пре- 
образование изображения из градаций серого в чёр- 


но-белое, расчёт ширины годичных колец и инте- 
грального прироста по всем азимутам, нормирова- 
ние ширины прироста каждого кольца. 

Процедура медианной фильтрации изображения 
использовалась для устранения импульсного шума. 
В результате медианной фильтрации происходит 
ограничение сигнала сверху и снизу отрезками, па- 
раллельными оси абсцисс. Это проиллюстрировано 
на рис. 3, где изображена функция плотности древе- 
сины до и после применения медианной фильтра- 
ции с пятиэлементным окном. Так как исходная 
функция плотности древесины, наряду с импульс- 
ным шумом, содержит и низкочастотный шум, то 
его устранение осуществляется фильтрацией в поли- 
номиальном базисе на основе метода наименьших 
квадратов. Из рис. 4 видно, что функция плотности 
древесины после применения полиномиальной 
фильтрации (полином пятой степени) не ухудшила 
своего вида, а низкочастотный шум исчез. В резуль- 
тате фильтрации получается сигнал, из которого вы- 
делена низкочастотная составляющая и импульсный 
шум при сохраненной информации о переходных 
зонах и незначительном изменении амплитуды. 



Рис. 3. Медианная фильтрация сигнала 



Рис. 4. Полиномиальная фильтрация сигнала 
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Также проводились вычисления значений ин- 
дексов прироста по ширине годичных колец дере- 
вьев - сканирование изображения по выбранному 
сечению, обработка полиномиальным фильтром, 
преобразование из градаций серого в чёрно-белое, 
подсчета средних значений по индексам прироста 
годичных колец, рис. 5. 

На графике по оси х отложены года, по оси у - 
прирост и значение индекса. Исходная ширина го- 
дичных колец деревьев отмечена квадратами в узло- 
вых точках, сглаживающая кривая - вертикальны- 
ми прямыми, индексы ширины - окружностями. 



Год 

Рис. 5. Вычисление значений индексов прироста годичных 
колец деревьев 


В данном случае индекс прироста - трансфор- 
мированное значение ширины годичных колец, по- 
лученное путём деления её на соответствующее ей 
значение на сглаженной кривой. Вычисление значе- 
ний индексов прироста осуществляется с целью 
устранения возрастного тренда из исходных измере- 

СПИСОК ЛИТЕРАТУРЫ 

1. Биоиндикация загрязнений наземных экосистем / Под ред. 
Р. Шуберта - М.: Мир, 1998. - 350 с. 

2. Розенберг Г. С. Биоиндикация: теория, методы приложения. - 
М.: Наука, 1983. - 163 с. 

3. Шиятов С. Г., Ваганов Е.А., Кирдянов А.В., Круглов В.Б., Мазе- 
па В.С., Наурзбаев М.М., Хантемиров Р.М. Методы дендрохро- 
нологии. Ч. I. Основы дендрохронологии. Сбор и получение 
древесно-кольцевой информации. - Красноярск, 2000. - 80 с. 

4. Ваганов Е.А., Шашкин А.В. Рост и структура годичных колец 
хвойных. - Новосибирск: Наука, 2000. - 232 с. 


ний путём перевода последовательности значений 
прироста в безразмерные значения индексов шири- 
ны годичных колец. В изменчивости индексов влия- 
ние внешних условий максимально и процедуры 
дальнейшего анализа позволяют выделить ведущие 
факторы среды и оценить их относительный вклад. 

В целом, полученные оценки позволяют произ- 
водить реконструкцию многих важных климатиче- 
ских характеристик, основными из которых явля- 
ются температура воздуха в различные сезоны и за 
год, количество осадков в различные сезоны и за 
год, аномалии атмосферного давления, повторя- 
емость и интенсивность засух, повторяемость за- 
морозков в течение вегетационного периода и 
сильных морозов в зимнее время. 

Заключение 

Представлена архитектура многомашинного 
комплекса, обеспечивающая динамическое изме- 
нение инфраструктуры системы обработки и до- 
ступ к программным средствам математического 
анализа годичных колец деревьев с целью монито- 
ринга окружающей среды. Математическая обра- 
ботка основывается на представлении плотности 
древесины как пространственно-временном коле- 
бательном процессе и включает вычисление азиму- 
та и среднеквадратичной ширины области макси- 
мального прироста, индексов прироста ширины 
годичных колец деревьев. 
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